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We study one dimensional wires with spin-orbit coupling. We show that in the presence of 
Zeeman field and strong electron-electron interaction a clean wire may form fractional helical liquid 
states with phenomenology similar to fractional quantum Hall liquids. Most notably, the wire's two 
terminal conductance is predicted to show fractional quantized conductance plateaus at low electron 
density. When the system is proximity-coupled to a superconductor, fractional Majorana bound 
states may be stabilized. We discuss how disorder destabilizes these fractional phases. Possible 
experimental realizations of similar states in double wire systems are discussed. 
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Introduction: Two-dimensional topological insulators 
have one dimensional helical modes at their edges. In 
these edges counter-propagating modes are related by 
time reversal [TJ [5] . These systems are most easily un- 
derstood when viewed as integer-quantum-spin-Hall sys- 
tems, in which electrons with opposite spins are at integer 
quantum Hall states with filling v = ±1. When properly 
subjected to proximity coupling to superconductors and 
ferromagnets these edges may host localized zero energy 
Majorana modes [3]. 

Similar phenomenology occurs in one-dimensional 
(ID) wires where electrons are subjected to spin- 
orbit coupling and an external Zeeman field [4]. The 
two-terminal conductance of such a wire varies non- 
monotonically as a function of increasing chemical po- 
tential, changing from zero (for a fully depleted wire), 
to a 2e 2 /h plateau, then to a e 2 /h plateau at which the 
wire is helical, and back to 2e 2 /h (see Fig. 1) as observed 
in a recent experiment [5], Furthermore, it was shown 
theoretically [H [3] and experimentally [3-^ that in prox- 
imity to a superconductor the helical state gives rise to 
end modes that are Majorana modes. 

In two-dimensions (2D) interactions between electrons 
may give rise to fractionalized phases, whose edges 
are gapless Luttinger liquids. In this work we study 
how this fractionalization may be reflected in quantum 
wires. Considering the same setting used for the one- 
dimensional helical quantum liquids in ^ E] , namely a 
single-mode wire where electrons are subjected to spin- 
orbit coupling and a Zeeman field, we examine the con- 
ditions under which one-dimensional "fractional helical 
liquids" may be formed, their properties and their sensi- 
tivity to disorder. The crucial dimensionless number for 
our analysis is the ratio of the wire's Fermi momentum 
k-p to the momentum characterizing spin-orbit coupling, 
kso- We find that when this ratio is an odd integer 2n+l 
a fractional helical liquid state may form (n is a non neg- 
ative integer). In such a liquid the two-terminal con- 
ductance of the wire is ^- i+pn+IF ' ^' ie non ~ mon °tonic 
dependence of the conductance on chemical potential be- 



comes then much richer. 
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FIG. 1: Left panel: spin-orbit split bands in a quantum wire, 
and a small energy gap opening due to magnetic field near 
the band-crossing point. Right panel: a sketch of the zero 
temperature conductance dl/dV versus gate voltage. Inset: 
two probe measurement setup where the fractional state is 
stabilized in the central region between two non-interacting 
leads; dashed lines represent the scattering pattern of bosons 
for n = 0. 



Furthermore, in 2D counter-propagating edge modes 
may occur in fractionalized phases as well. Examples 
range from fractional topological insulators, that are 
presently theoretical, to graphene electron-hole bi-layers, 
which are experimentally feasible. It was shown theo- 
retically that in proximity to superconductors and fer- 
romagnets the counter-propagating edge modes may be 
gapped to host localized fractionalized Majorana modes, 
that are non-abelian anyons of a new type 10 13J . We 
show here that proximity-coupling of a ID fractional he- 
lical liquid to a superconductor gaps its spectrum, and 
gives rise to fractionalized Majorana modes akin to those 
found in 10 \'A . The resulting gapped state is not within 
the list of possible gapped one-dimensional phases clas- 
sified in Refs. Q31 [T5] It then comes as no surprise that 
these states are not stable to disorder. 

Our analysis benefits greatly from the study of the 
two-dimensional quantum Hall effect in terms of a set 



2 



of one-dimensional coupled quantum wires. In this ap- 
proach, carried out in [16] and [17], the 2D plane is seen 
as made of an array of coupled wires. The spectrum of 
each wire is parabolic, and the parabolas are shifted with 
respect to one another in momentum space by the mag- 
netic field. The number of wires is proportional to the 
sample's width, and approaches infinity in the thermody- 
namic limit. As we show below, our system corresponds 
to the system considered in Ref . Lo] when the number of 
wires is two. 

While the main focus in this manuscript is on a single 
wire with spin-orbit coupling that shifts the free electron 
spectrum parabolas, a similar effect may occur in a pair of 
parallel wires subjected to a magnetic field perpendicular 
to their plane. In these wires, spin polarization by the 
magnetic field and a proper tuning of the Fermi level in 
each wire may stabilize fractional helical phases. Such 
wires were experimentally realized by means of cleaved 
edge overgrowth [T8] . 

Model: Our model of interest consists of a spinfull, sin- 
gle band wire, parallel to the x direction, with Rashba 
spin-orbit (SO) coupling u, and with a magnetic field B 
nonparellel to the SO quantization axis. The Hamilto- 
nian is 
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with a\ i = x,y, z being the Pauli matrices (here and 
henceforth we work with a system of units in which 
h = 1), and k being the momentum. The SO coupling 
causes a shift of the momenta k of the free electrons 
parabolic spectrum E(k). The spin up electronic dis- 
persion (up with respect to a direction y set by the SO 
coupling) is shifted to the left by a momentum k^o = mu 
while spin down is shifted to the right by fcso (see Fig. 1). 
The two parabolas intersect at zero momentum so that 
an application of uniform (Zeeman) magnetic field (not 
parallel to y) opens a gap in the spectrum near zero mo- 
mentum. Two electron helical modes, one propagating 
to the left and the other to the right, remain free. 

Away from the gap there are four Fermi points [with 
wave vector ki = ±(fcgo ± &f)] and therefore we define 
4 fermionic fields ipi with I = 1,2 corresponding to the 
left and right moving modes of <r y = 1 and I = 3, 4 corre- 
sponding to the left and right moving modes of o~ y = —1 
(see Fig. [I). 

Interactions affect this Hamiltonian in two ways. Small 
momentum interactions are accounted for by introducing 
Luttinger liquid parameters, as usual. Large momentum 
interactions need to be elaborated on for our discussion. 
As pointed out by Kane et al. in Ref. [16] in the context of 
the fractional quantum Hall (FQH) effect, multi-electron 
processes that involve large momentum transfer may lead 
to an opening of an energy gap at proper values of the 
chemical potential. Following [16] we focus on the inter- 
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For n = this term reduces to the mixing of modes 2 and 
3 at the band crossing point, Oq — gs^ipl- For n = 1 
this process takes 3 particles across the 4 Fermi points, 
2 — > 1, 2 — > 3, 4 — > 3; see arrows in Fig. ([I]). Such a term 
(together with its hermitian conjugate) is odd under time 
reversal, hence it requires magnetic field to be present. 
It is generated at second order in the bare interaction 
strength gg oc -Bf/| fe , with U q the interaction potential. 

The combined Hamiltonian is most conveniently 
treated by bosonization. The four chiral bosonic fields 
corresponding to the four fermionic fields are tpi — 
e 1 ^ 1 ; tpjtpi — (— l) l j^d x (f>i. The bosonic Lagrangian 
density is 



+ g B COS {03 - 02 + n [(03 - 04 ) 



^2-0i)] + fca;} 
(3) 



The positive definite symmetric matrix U, describing the 
bare velocities of the fields and the low momentum trans- 
fer interaction, is not universal, here k = k^ — k2+n[(Jc3 — 
ki) — (k 2 ~ ki)]. When k « 0, the term proportional to 
<7b gaps out two of the four modes in the spectrum. Us- 
ing the parabolic spectrum when B < Ajo = mu 2 /2 
we find that k = for k F = k so /(2n + 1). With 
the relation p = 2kp/ir, where p is the density of the 
spinfull wire, we find that a gap will open for filling 
v = Ttso ~ ~ 2^-4-1 • Notice that a component of 



2n+l ' 



the Zeeman field in the SO coupling direction (y) will 
change the conditions for k = 0, the generalization is 
straightforward. 

The gs term in the Lagrangian ^ motivates us to 
introduce a new set of chiral bosonic fields, denoted by 
m = Ysj Mijfa, S iven explicitly by: 
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Notice that for n = 0we have r\i — <pi . In terms of the 
r\ fields, the Lagrangian density becomes, 

i ij 



fjB 



cos(r,3- ??2 )--i^ 2 ^ TI??i , (5) 
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with U = {M- l ) T UM- 1 . 

For B 3> Ago the description in terms of four bosonic 
fields ceases to be valid (as k 2 and k 3 do not exist any- 
more). However, similar to the case for n = when the 
system is fully gapped (e.g. by proximity to a supercon- 
ductor) the gap is not necessarily closed when B > Aso 
and a phase transition would not necessarily occur then. 

We now turn to discuss several properties of the wire 
that follow from the analysis of Lagrangian |5]) in the 
limit of large gs- 

Two terminal conductance: In a two terminal con- 
ductance measurement one connects two noninteracting 
reservoirs to the wire. We model this situation by a cen- 
tral interacting section connected to the reservoirs via 
two non interacting sections. By assumption, the transi- 
tion between the sections is smooth on the scale of kp so 
that a 2/cf scattering does not occur. The DC conduc- 
tance can then be found by studying a scattering prob- 
lem in which electron currents incident from the reser- 
voirs scatter in the central region as shown in the inset 

to Fig.rj 

To solve the scattering problem it is useful to define 
incoming and outgoing current vectors. The incoming 
currents are T = (Ji, I 2 , 13, h) T , with I x = j^d t $i\ x=Xl 
and I = 1, ... ,4. The points xi are at +00 for the left 
movers I = 1, 3 and — 00 for the right movers I — 2, 4 (see 
Fig. A similar definition holds for the outgoing cur- 
rent O = (Oi,0 2 ,0 3 ,0 4 ) , with xi — s- —xi. Assuming 
now that incoming left movers (1,3) both emanate from a 
reservoir at potential V and the right movers (2,4) from a 
reservoir at potential zero, we have I = \V (1, 0, 1, 0) T , 
and O = SI, with S being the scattering matrix. The 
conductance G is then given by: 

GV = h+h-0 2 -0 A = -h-h + O z + O l 
G/(e 2 /h) = 2 -(0,1,0, 1)5(1, 0,1,0) T 

= (1,0, 1,0)5(1,0, 1,0) T . (6) 

Next we should find the scattering matrix S. In the free 
case gs — and the scattering matrix is equal to the 4x4 
identity matrix. Simple algebra then gives: G = 2e 2 /h. 
For the n — case with gg — > 00 the chiral field fa is 
fully reflected to fa (see inset to Fig. [I]), and this leads 
to a scattering matrix Sq equal to a matrix in which the 
second and the third rows of the rank 4 identity matrix 
are interchanged. For that case G — e 2 /h as expected. 

Determination of the scattering matrix for n =/= 
follows directly from the structure of the Lagrangian 
Eq. ([5]). For gs — > 00 the cosine term implies that 
V2( x ) = V3( x ) f° r ~~ L/2 < x < L/2 which is the region of 
the wire where the electrons interact. Using this equation 
at x = —L/2 and at x = +L/2 gives two relations 



(n + 1)Z 2 - nOi 
(n + l)0 2 - nil 



(n + l)0 3 - nh, 
(n + 1)J 3 - nO A . 



A second pair of equations is obtained by realizing that 
fields 771 and 7/4 propagate freely through the interacting 
wire implying r]x ! i(—L/2) = T]± t i(L/2). In terms of the 
original fields the equations for r]i and 774 become 



(n + l)0 1 -n0 2 = (n + l)h-n0 2 , 
(n+l)h-n0 3 = (n + l)0 4 -nl 3 , 



(8) 



respectively. This set of four equations can be written in 
a matrix form as AO = BI where A, B are 4x4 matrices, 
giving the desired result S = A~ X B. Substituting this 
result for S in ^ we find the two terminal conductance: 
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h 1 + (2n + l) 2 



hv 2 + l 
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For n — we find G — %- as expected, for n = 1 
1 2 

[y = 1/3) G = g ^- and so on. In general, unlike the FQH 
effect the conductance is not equal to the filling factor v. 
Similar to the two modes on the edges of a FQH 2D sys- 
tem, in the present case two modes remain untouched by 
the massive term that mixes 7/ 2 and 773. However, unlike 
the 2D system were the free modes on different edges 
emanate from different reservoirs and propagate freely, 
in the present case the free propagating modes in the in- 
teracting section are mixed at the boundary between the 
interacting and the noninteracting sections. The physical 
proximity of the two counter-propagating modes makes 
them susceptible to disorder, as we will discuss below. 
We note that as in the case of a Luttinger liquid describ- 
ing quantum Hall effect or quantum wires, the low mo- 
mentum interactions Uij do not affect the two terminal 
conductance [TM2"3"] . 

Refs. [ini [T7]show that in the limit of an infinite num- 
ber of spinless wires whose parabolic dispersion is shifted 
by a magnetic field the two-terminal conductance ap- 
proaches the expected value of the quantized Hall con- 
ductivity. The case we study corresponds to two wires. In 
the appendix we study numerically the case of N wires, 
and show that the deviation of the two terminal con- 
ductance from the quantized Hall conductivity decreases 
very quickly with N, and is the order of e~ N for N > 5. 
This reflects the fact that the tunneling between the edge 
modes is exponentially small. All the values that we find 
for the two terminal conductance of N wires at filling 
v = kp/kgo = l/(2n + 1) conform to the expression 
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Measurements of individual elements of the S -matrix: 
The individual elements of the 5-matrix may be mea- 
sured in a setup such as used by Auslaender et. al [18] . In 
this setup two parallel wires are subjected to a magnetic 
field perpendicular to their plane. If the spins are po- 
larized and tunneling between the wires is enabled, this 
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system maps onto the system we study here, with the 
wire index playing a role of a pseudospin. The two wires 
are connected to four separate reservoirs, allowing a sepa- 
rate control of the incoming currents 1% . . . J4. For exam- 
ple, the application of a voltage difference V between the 
wires corresponds to an incoming current / = (V, V, 0, 0), 
and results in a net current flowing between the wires. 
This current is I s = G S V = (0, 0, 1, 1)S(V, V, 0, 0). Using 
the above S matrix this current is found to be equal to 
the charge current that flows when voltage is applied be- 
tween the right and left moving modes, i.e., G s = G. In- 
terestingly, even when the current from wire 1 emanates 
equally from the two reservoirs the wire is connected to, 
the magnetic field between the wires results in a cur- 
rent partitioning in wire 2 that is not symmetric between 
the left and right drains. Rather, the difference between 
the currents flowing to the left and right in wire 2 is 



(0,0, 1,-1)S(V,V, 0,0) 



_ e 2 2(2n+l) 
h l+(2ra+l) 



V. 



Disorder: Disorder that leads to scattering between 
the free left and right movers spoils the quantized con- 
ductance. This is obvious in the n = case, where in- 
teractions are absent. For n 7^ single electron tun- 
neling is of the form tpjipi'. In the bosonized form it 
is A;/' cos(0i — 4>i>). To find their relevance we write 
them in terms of the 77's. Since 772 is pinned to 773 
by the gs term, 772 + 773 fluctuates wildly. As a con- 
sequence, a term that contains 772 + 773 is irrelevant, 
while the combination 772 — 773 can be approximated as 
a constant. These considerations make the A14 and A23 
the only relevant single particle terms. These terms 
translate to cos [a(?7i — 774) + /3(?72 ~ ^3)]) where (a, j3) — 
2^+i(n+ 1,77) for A 14 and (a, (3) = ^+i( n ' n + x ) for 
A23. Both terms scatter charge between the two com- 
pressible modes 771 and 774, and affect the two terminal 
conductance. 

Generally, the scaling dimensions for these two terms 
makes them relevant. For repulsive interactions, these 
dimensions are bounded from above by 2a 2 /(l + v 2 ), see 



Eq. (32) in the appendix. However, the divergences as- 



sociated with the relevance of impurities are cut-off here 
by the finiteness of the wire, as well as by finite voltage 
and temperature. 

Multi-particle scattering by impurities may be even 
more relevant than single particle scattering. For exam- 
ple, Akp scattering of the form cos (</>i — + (f>3 — 4>i) = 
cos ^(771 — 774 — 772 + 773) leads to a scaling dimension 



bounded by 2t/ 2 /(1 



This backscattering opera- 



tor corresponds to a quasiparticle tunneling in the FQH 
2D limit. 

Fractional charges: In the n = case a soliton (i.e. 
a kink between 0, 2tt in the 772 — 773 field) describes an 
electron or hole with charge C = ±e. Since l)'<^z — 

2w 1 +1 1) vi it i s readily understood that for general 



71 a kink in the the 77 fields carries a charge [Trj] 

e 



C = ± 



277 +1 



We note, however, that in ID those quantized charges 
can not be isolated since the free modes (composed of 771 
and 774) will screen those charges in a way that depends 
on the low momentum interaction Uij. 

Fractional Majorana modes: So far we considered a 
compressible helical liquid, in which the fractionalized 
nature of the wire reflects itself in the fractional two ter- 
minal conductance. The interaction term gs cos(?73 — 772) 
gaps two of the four gapless modes in the wire, but leaves 
the other two gapless. By coupling the wire to a super- 
conductor through the proximity effect, the wire may 
be fully gapped. In the 77 = case a coupling to su- 
perconductor leads to the formation of Majorana end 
modes [H [5] . We will show now that in the n > 1 case 
proximity to a superconductor leads to the formation of 
fractional Majorana modes. 

The proximity induced superconducting mean field 
Hamiltonian projected to the four Fermi points is "Ha = 
A J dx(ipi(x)ip4 : (x) + tp2(x)ip3(x) + H.c.). Bosonizing, we 
find Ha ~ Acos(0i + 4>\) + Acos(</>2 + <fe)- Both op- 
erators are irrelevant for n > 1 since they contain an 
exponent of 772 + 773 which is conjugate to a pinned field 
and hence its correlation function decays exponentially. 
We thus look for an operator that couples to the super- 
conductor but does not contain 772 + 773. The leading 
operator of this form is 

c£ = gA^M^n^tr^l ~ 9a cos(7 ?1 + 774). (f 1) 

(In the wire the operator cx cos (771 — 774) does not con- 
serve momentum and hence is absent in a translationally 
invariant system.) This process is shown in Fig. J2jj. For 
77 = I (at filling factor v = 1/3) it is generated at first 
order in A and in the interaction, g^ oc A[/| fe . Thus 
we see that when the proximity coupling is strong enough 
this term will pin 771+774 and hence the low energy theory 
will become gapped. A similar many body process will 
also generate a term cos(?72 + 773) which competes with 
the gs term, see Fig. [2^,. As long as the latter is stronger, 
however, we may neglect this term. 

We now claim that the resulting gapped phase has frac- 
tional Majorana modes at the ends of the wire. To treat 
the boundaries of the finite length (L) wire we use the 
usual unfolding transformation [24s which results in a 
boundary condition 

0i(as) = Mx), Mx) = M*)> ( x = 0,i). (12) 
In the new variables the boundary conditions become 

m(x) = V2(x), 773(2;) = 774(2;), (x = 0,L). (13) 
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FIG. 2: (a) Superconducting analog of the process in Fig. 1. 
(b) Many particle process inducing a superconductivity gap 
in the fractional helical liquid. 



It proves convenient to define a theory with a single 
left moving field £l(x) and a single right moving field 
£r(x) on a ring of circumference 2L, according to 



£r(£) 



(x)= U(3)(*), 



< x < L, 



1 ?7i(4)(2£-x), L<x<2L, 
see Fig. ([3]). The magnetic and superconducting pertur- 








FIG. 3: Unfolding transformation on the spinfull finite wire 
to a ring whose upper part is coupled to the magnetic pertur- 
bation and whose lower part is coupled to the superconducting 
perturbation. Fractional Majorana modes form at x = and 
x = L. 



bations to the Hamiltonians now act in different regions, 



5C = 



2L 



dxg B (x) cos(f fl - £ L ) + g A (x) coa(£ R + dJ;L5) 



where gs is finite in the region < x < L and g& is 
finite in the region L < x < 2L. For n — the interface 
between two gapped regions dominated by cos(£# — £l) 
and cos(£a + binds Majorana zero modes. For n ^ 
0, the interfaces bind fractional Majorana modes as was 
recently demonstrated in [TUHT3"] . 

In principle, such zero modes lead to a degeneracy of 
the ground state. However, this degeneracy is apparent 
only when the system hosts at least four zero modes, 
such that there are several ground states that share the 



same fermion parity. Systems that host four zero modes 
may be double wire systems, or a single wire divided to 
segments of phases that are topologically different. In 
the present context, we may imagine constructing such 
segments by having regions where the <7b,<?a are opera- 
tive, alternating between regions where c/a is put to zero, 
and the r}i , 774 modes are gapped by a perturbation of the 
form 



9B> 



cos (rji — 774) dx 



(16) 



where a, b are the end-points of the segment. Physically 
this term corresponds to a spatially oscillating magnetic 
field with wave vector ikso- Alternatively, four zero 
modes may be constructed in systems of two wires with 
properly-tuned couplings that allow also for a transfer of 
charge between the wires. 

Disorder has a destabilizing influence on the 
gapped superconducting phase and its degener- 
acy. Indeed, a backscattering term of the form 
cos [a{i]i — 774) + /3(?72 — 773)] will lift the degeneracy, 
as now for example the states with 773 — 7/2 =0 and 
773 — ^72 — 27r have different energies. 

Relevance and irrelevance: The fractional helical 
states, as well as the superconducting phases, can be sta- 
bilized if both gs and c/a are large enough. But yet a 
relevant question is under which conditions it is possible 
to trigger those instabilities even when the bare pertur- 
bations gs and g& are small. We will use the renormal- 
ization group (RG) method to address this question. We 
note, however, that for finite one-dimensional wires, the 
RG flow is cut-off by the wire's length, and does not nec- 
essarily get to its asymptotic value. Therefore, the rel- 
evance of perturbations to actual experimental systems 
will strongly depend on the bare values of couplings. 

In an interacting Luttinger liquid described by charge 
and spin interaction parameters K c and K s , respectively, 
we find that the operators [in Eq. ^] and [in 
Eq. (11)] have scaling dimension 



x B = \{K~ X + (2n + l) 2 K c ), 
x A = \{K~ X + (2n + l) 2 K s ), 



(17) 



respectively. For small magnetic field we have K s = 1 
due to an SU(2) symmetry of model ([!]) realized in terms 
of fermions ip'i(x) = e lStxkso ipi(x) with si = s-2 = —S3 = 
-a 4 = 1 [25]. 

Consider first the magnetic perturbation gs- If the 
density is fixed at v = l/(2n + 1), then the instabil- 
ity into the fractional helical state occurs for infinitesi- 
mal magnetic field when the scaling dimension is equal 
or smaller than the space-time dimension, namely at 
Kc < 3/(2n + l) 2 . If fcp deviates from the special value 
fcso / (2n+ 1) then a non zero bare value of gs is essential 
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to stabilize the helical phase through a commensurate- 
incommensurate (C-IC) transition. As known from the 
study of C-IC transitions [24] at the transition the di- 
mension xb takes the value one. Repulsive interactions 
reduce the value of the critical (bare) qb [2U] . 

A rough estimate of K c in wires shows that it may 
actually be rather small. To linear order in the inter- 
action U q the charge Luttinger parameter is given by 

(4C7 \ — 
1 + vh J J . Since the Coulomb interaction 

U q diverges logarithmically with the distance to the gates 
and near the band bottom the density of states 1/vf di- 
verges as well, their combination leads to a small K c . 

The superconducting perturbation g& does not oscil- 
late in space and hence does not require a special density. 
However, except for n = 0, it is irrelevant, xa > 2. In- 
deed, even in the n = case it was pointed out that the 
superconducting gap is reduced due to interactions [27] , 
Thus, strong proximity coupling is necessary. 

We note that apart from the couplings gs and (/a, the 
system may become unstable with respect to spin density 
wave formation Osdw — V 7 ! ■ Note that there is 
no special filling condition for this operator to conserve 
momentum and this operator is exactly marginal for all 
values of K c . The strong coupling picture of this instabil- 
ity corresponds to an antiferromagnetic spin arrangement 
with wave vector 2k p. This instability can be inhibited 
by applying a magnetic field parallel to the y axis, shift- 
ing the up-spin and down-spin bands vertically. When 
that happens the SDW term does not conserve momen- 
tum since at the chemical potential the spin up Fermi 
momentum fep-f = (&2 — ki)/2 is not equal to the spin 
down Fermi momentum fcp^ — (^4 ~ k\)j2 and the cor- 
responding SDW phase can not be stabilized. 

Summary: In this manuscript we showed theoretically 
that clean spinfull wires with spin-orbit coupling and 
strong electron-electron interaction may give rise to novel 
fractional helical liquids. The two terminal conductance 
of these liquids is quantized to a fractional value. We 
argued that in proximity to a superconductor fractional 
Majorana modes may be stabilized. We discussed the 
influence of disorder on these wires and pointed out that 
alternative systems such as the two coupled wires (that 
do not have substantial spin-orbit effects) may give rise 
to a similar fractionalization. 

We expect that the ideas put forward here for fraction- 
alized one-dimensional systems coupled to a supercon- 
ductor can be extended, by studying several wires con- 
nected in parallel, to novel two-dimensional fractionalized 
superconducting phases. Furthermore, it was shown in 
Ref. [T7] that two dimensional non-abelian quantum Hall 
systems whose quasi-particles allow for universal topolog- 
ical quantum computation may be viewed as constructed 
out of an array of parallel wires, with properly tuned 
inter-wire tunneling and interactions. We envision that a 
combination of such interaction, together with spin-orbit 



coupling and proximity coupling to a super-conductor 
may allow for such computation to be carried out in clean 
quasi-lD systems as well made of several wires. 

We acknowledge discussions with J. Alicea, O. Aus- 
laender, M. Barkeshli, E. Berg, D. Clarke, L. Fidkowski, 
M. Freedman, N. Lindner, and A. Yacoby. The work was 
supported by a BSF, ISF, DFG, Minerva and Microsoft 
grants. 



APPENDIX 
Generalization to TV parallel wires 

Kane et. al. [TB] studied the 2D FQH system as many 
parallel spinless wires where a term of the form Eq. |2]) 
couples each pair of neighboring wires. We will now focus 
on such a system with a finite number of wires, N. This 
analysis will bridge between our system of interest, of a 
spinfull wire with SO coupling (N — 2) and the 2D limit 
(N oo). 

Each wire j — 1,...,N is described by left (cf>2j-i) and 
right (4>2j) boson fields. The chiral fermion operators are 
given by %l)j(x) oc e^ [x) . The La grangian under consid- 
eration is C = Cq + 5C, where we assume that E7y of 
Eq. ([5| is diagonal and 

- 2N 

£ o = ^ 12 ((-iWA^ + v F {d x ^f + 2E<i> j (-iy) . 



By a straightforward generalization of Eq. Q we define 
r/j fields, j = 1, 27V, in terms of which the perturbation 
5£ reads 

N-l 

6C=^2g B cos (?72j - V2j+i) ■ (19) 
i=i 

The fields rji and 772^ will form the low energy sector 
providing the edge theory. 



Two-terminal conductance 

We now consider the two-terminal conductance in the 
case that the interactions act in a finite region —L/2 < 
x < L/2. The analysis below is a direct generalizations 
of Eqs. ([7]), and |8]). Consider a situation where 
the non-interacting lead to the right has chemical poten- 
tial V, and the lead to the left is at chemical potential 
0. We now describe the scattering problem in terms of 
incoming fields 1 L (i = 1,...,N) and outgoing fields L 



7 



(i = 1, N). The conductance is given by 



N 



G N V = 53(/a-i - O w ) 



z=i 



W-(0,1,0,1,...)S(V, 0,V,0,...) 



= 5Z(-/2; + 2; _!) 



z=i 



(1,0,1,0,...)S(V,0,F,0,...). 



(20) 



The 27V x 27V scattering matrix O = SI needs to be 
determined by a set of 27V equations. Generalizing the 
procedure in the main text, we get (27V — 2) equations 
by setting r)2i — f]2i+i at x = — L/2 and at x = L/2, 

(n + - n0 2l -i = (n + l)0 2 i+i - nI 2i+2 , 

(n + l)0 2 i-nl 2 i-i = (n + l)I 2l+1 - n0 2l +2, (21) 

where i = 1,2, ...,N — 1. Two additional equations orig- 
inate from the free propagation of r\\ and r\2N 

(n + l)h-n0 2 = (n+l)Ox -nl 2 , 
(n + 1)I 2N - n0 2N - X = (n + l)0 2N -nI 2N -i.(22) 

This set of equations can be written in matrix form AO = 
BI , giving S = A^ 1 B, and the two terminal conductance 



is found using Eq. ( 20 ) 



The result of this calculation is that the conductance 

Explicitly, we find that 



(2n+l 



exponentially tends to 
the result is given by Eq. ( 10 ) for all v and TV, which is 
plotted in Fig. g| for n= 1 (u = 1/3). 




FIG. 4: Two terminal conductance on an array of TV wires 
for n = l(u= 1/3). 



We remark that the sum of each column of the S 
matrix is equal to unity, implying current conservation. 
However, S can have negative elements. For example, for 
TV = 2, 



S = 



1 \ 



2 14-2 
-2412 
\ 1 -2 2 4 j 



(23) 



For an incoming pulse (soliton) at I\, we have as out 
going currents O = |(4, 2,— 2,1). The negative current 
in 3 is not forbidden as it corresponds to an antisoliton. 



Low energy action 

We now derive the low energy effective action. This 
derivation will lead explicitly to the two freely propa- 
gating chiral modes as specific linear combinations of r\\ 



and r\2N , see Eq. ( 28 1 . For finite TV this effective theory 



is different from the edge theory of a FQH state, as dic- 
tated by the parameter Kjq derived below, determining 
the commutation relations of the chiral fields. 

In terms of the rj fields the free Lagrangian becomes 



2N 

j=i 
vpC 



Ev 



8tt 



4tt 

2 VfD 



4?r 



2tt 

d x r) 2j -id x r) 2j ] 



(24) 



with C(D) = 1 ± v 2 . We see that for v ^ 1 there is 
coupling between fields f] 2 j-i and 772^ that will medi- 
ate a coupling between 771 and 7727V • The effective action 
for qb — > 00 is obtained after integrating out the fields 
{{V2j + V2j+i), 3 = 1, TV - 1} and setting rj 2j = rj 2j+1 
(j = 1,...,N-1). Writing Eq. Q in terms of rfcj+Tfcj+i, 
and setting r\ 2 j — 772.7+1 = we find 



j=l,j=2JV 

vpC 



An 



-d t rijd x T)j 



Ev 
2^ 



8tt 

vf 
16tt 

v F D 

~8tT 



(dxVjY 



N-l 



£ 9 x(V2l + V2l+l)Uljd x (f] 2 j + V2j+l) 
1,3=1 

(d x (r]2 + 773)^771 + d x (r\ 2N - 2 + r] 2 N-i)d x T} 2N ) , 



where the (TV — 1) x (N — 1) matrix U is given by 

\ 



U 



(CD/2 

D/2 C D/2 

D/2 C 

V = 



(25) 



/ 



To get the interaction between 771 and rj 2 N we complete to 
a square the terms containing d x (772+ 773) and d x {i]2N-2 + 
V2N-1) and obtain the effective Lagrangian 



£ e ffD?l,772-/v] 

Ev 



V 



4tt 

v f 



(d t i]2Nd x r] 2N - d t 7]id x 7]\) 



., (r? 2 jv - m) + t£- [(C + FiWxim + m N )) 2 
2tt 1o7t 



(C + F 2 )(d x ( Vl - rj 2N )f 



(26) 



with 



Fi, 



(27) 



JV-1,JV-1 



Finally, to bring the action to a non-interacting form, we 
define new propagating eigenmodes 



2 " L 'l/Vv ") " JV (a'v 



and obtain 



1 

4tt 



(28) 



(29) 



Here p = R/L = ±1. i^T/v and u are given by 



1 

Kn 



F 2 



= nu- 1 y /(C + F 1 )(C + F 2 ).(30) 



Here u is identified as the velocity. For N = 2 the ma- 
trix U has one element, C, and one obtains K 2 — v - 2 +1 ■ 
Generally we find that this definition of Km coincides 
with the result for the two terminal conductance K N l — 
Gn/{s 2 IK) [where the general equation for Gat is given 
by Eq. (10)]. Thus, also K^ 1 tends exponentially to the 
value v, where the effective theory reduces to two inde- 
pendent chiral FQH edge theories [28] . The total density 
* s P = ^^x4>R ~ ^x4>l- I n the DC limit variation with 
respect to d x <fi p gives the compressibility is -£ — ■ 
Perturbations to the helical liquid phase of the quan- 
tum wire, such as impurity scattering, can be studied via 
renormalization group analysis with respect to the fixed 
point of Eq. p0| (and setting N = 2 for the case of a 
quantum wire) . The scaling dimension of the qua sipa rti- 
cle operators e l<f>p is (p = L/R). Using Eq. (28) we 
see that 



v(j]\ - r) 2N ) = 4>l~ 4>r- 
Hence the dimension of cos [a{rj\ — 772 jv)] is 



(31) 



(32) 
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